Equilibrium geometries of low-lying isomers of some Li clusters, within Hartree-Fock 
theory plus bond order or MP2 correlation corrections 
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In a recent study by Kornath et al. [J. Chem. Phys. 118, 6957 (2003)], the Li n clusters 
with n — 2,4 and 8 have been isolated in argon matrices at 15 K and characterized by Raman 
spectroscopy. This has prompted us to carry out a theoretical study on such clusters up to n = 10, 
using Hartree-Fock theory, plus low-order M0ller-Plesset perturbation corrections. To check against 
the above study of Kornath et al., as a by-product we have made the same approximations for 
n = 6 and 8 as we have for n = 10. This has led us to emphasize trends with n through the Li n 
clusters for (i) ground-state energy, (ii) HOMO-LUMO energy gap, (iii) dissociation energy, and (iv) 
Hartree-Fock eigenvalue sum. The role of electron correlation in distinguishing between low-lying 
isomers is plainly crucial, and will need a combination of experiment and theory to obtain decisive 
results such as that of Kornath et al. for Lis. In particular, it is shown that Hartree-Fock theory 
plus bond order correlations does account for the experimentally observed symmetry T^ symmetry 
for Lis. 



I. BACKGROUND AND OUTLINE 

Lithium clusters by now have been studied by a variety 
of experimental techniques which include electron spin 
resonance, laser induced fluorescence, depletion spec- 
troscopy, photoionization, and Raman spectroscopy 1 . 
This Ref. [lj has investigated, in particular, the lithium 
clusters Li„, for n — 2,4 and 8. This has been done 
by isolating these clusters in argon matrices at 15 K. 
Whereas Kornath et o/i point out that most of the tech- 
niques listed above are difficult to correlate with the clus- 
ter geometries, their experimental work was able to de- 
termine the geometry of Li4 as a rhombic structure (O211) 
and for Lis a hypertetrahedral structure (T^) was shown 
to be in agreement with their Raman studies. 

This experimental work has motivated the present the- 
oretical study of Li„ clusters. Most attention is focussed 
here on n = 10, but we have also included n = 6 and 
11 = 8. In the latter case, we must mention the previ- 
ous theoretical work of Bonacic-Koutecky, Fantucci and 
KouteckjiS. 

The outline of the paper is then as follows. In Sec- 
tion [H] some general trends of the electronic structure of 
Li n clusters are depicted with n ranging from 1 to 10. 
We note that Kornath et al~ studied n = 2, 4, and 8 
extensively using high quality techniques. For Lig our 
more modest techniques essentially confirmed the find- 
ings of Kornath et al^ for the equilibrium geometry and 
the vibrational frequencies. So we focus in Section IIIII 
dominantly on Lig. Section llVI presents results for Liio- 
Finally, discussion of some additional trends, plus some 



suggestions for future work are covered in Section 



II. QUANTUM CHEMICALLY PREDICTED 
TRENDS IN Li„ CLUSTERS FOR n = 1 TO 10 



In this Section we present our results of the calcula- 
tions made on Li n clusters {n = 2, 4, 6, 8, 10), by using the 
Gaussian package (G03 Linux version)^. All the calcula- 
tions were performed using the standard self-consistent 
field Hartree-Fock (HF) theory with the 6-311G* basis 
set. The correlation energy was obtained using low-order 
M0ller-Plesset (MP2) corrections^, considering all (va- 
lence and core) electrons. For both HF and HF+MP2 
level, full geometry optimization and vibrational analy- 
sis were performed for all Li clusters. 

While Sections IIIII and Hvl will develop further the sys- 
tematic studies we have made especially for n — 6 and 
10, with the approximate quantum chemical approach set 
out above, this Section will present the general trends 
predicted by HF+MP2 theory for n = 1 to 10. 

To this end, Fig. ^ depicts the ground-state energy 
per atom for increasing n from 1 to 10. No commen- 
tary is needed on the results for n — 1 — 4, except to 
show as an inset the accepted planar structure of Li4. 
The HF energy of this cluster is —29.759 and the lower- 
ing by second-order M0ller-Plesset perturbation theory 
yields —29.859, both in Hartrees (as all energies below, 
unless stated otherwise). 

As seen in Fig. ^ t w0 structures were found to be 
low-lying isomers of Lig, with symmetries and D3^, 
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FIG. 1: Ground-state energy per atom of Li„ clusters, with 
n = 1 — 10. Inset: accepted D2h planar structure of L14. 
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FIG. 3: Low- lying isomers of Li6. Left: Planar structure, 
with F>3h symmetry. Right: Most stable isomer predicted by 
HF+MP2, with D,^ symmetry, and a dimer-like Li2 structure 
perpendicular to the plane. 



to consider our detailed results for Lig and, quite briefly, 
also Lig in Section ILUl immediately below. 



and these will be discussed in detail in Section lHfl below. 
Likewise, for Lig, C s and structures emerge and again 
Section lllll presents details. Li 10 has the lowest MP2 
energy per atom of the clusters considered and is the 
subject of Section ITVl 

Again concerning general trends, Fig.Ejshows the vari- 
ation of the HOMO-LUMO gap for Li„ clusters with the 
number of atoms n in the lithium cluster. It is of in- 
terest to note that though the variation is irregular for 
n = 4 — 8, the smallest gap is at n = 10. It is relevant 
to add here that solid lithium is metallic, with therefore 
zero energy gap. Of course the gap at n = 10 is still 
~ 5 eV! It is the trend, not the absolute numbers, which 
should be the focus here. 

With this relatively brief discussion of trends we turn 
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FIG. 2: LUMO-HOMO gap for Li„ clusters as a function of 
n. 



III. GEOMETRY, ENERGY, AND 
VIBRATIONAL FREQUENCIES OF Li 6 AND Li 8 
USING HF+MP2 THEORY 

In Ref . H 

a comment was made in their experimental 
considerations of a result which 'may indirectly support 
the absence of a Lig species.' Motivated by this, we have, 
by geometry optimization, 'converged' on the geometries 
of two low-lying isomers which are depicted in Fig. [3] 
The first of these is the planar structure (Fig. |3[ left) , 
with HF energy —44.658 and with addition of MP2 is 
—44.810. The present HF+MP2 approach, however, pre- 
dicts as the isomer we identify as the ground state, the 
planar square with a 'diatomic' Li2-like species perpen- 
dicular to the plane, through the center of gravity. While 
this has a (very) slightly higher HF energy than the pla- 
nar structure, it lies below that structure energetically 
when 'correlation' treated at the MP2 level is added, the 
results being E^if = —44.655 and with the MP2 correc- 
tion -44.824. 

Turning quite briefly to Lig, the geometry was settled 
as Trf from the Raman studies of Ref. [J whereas as can 
be seen from Fig. ^ the C s symmetry is predicted as 
very slightly lower by our present approximation. This is 
enough to illustrate the stringent test of many-electron 
approximate theories that is afforded by attempts to pre- 
dict decisively the ground-state energies of isomers of Li n 
clusters. However, the vibrational frequencies of Lig (T<j) 
are more in line with the observed values than those of 
the C s structure, which is somewhat encouraging in the 
light of the known geometry (T^) of the lowest lying iso- 
mer. 

We turn next to discuss, with fuller details, the Liio 
cluster. 
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FIG. 4: Low-lying isomers of Liio. Left: Asymmetric isomer. 
Right: Symmetric (D2d) isomer. The T>2d isomer of Liio is 
characterized by reflection symmetry with respect to the plane 
indicated in the figure, and by 180° rotation symmetry with 
respect to the C2 axis. See also- for calculated properties of 
both isomers. 
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IV. GEOMETRY, ENERGY, AND 
VIBRATIONAL FREQUENCIES OF TWO 
LOW-LYING ISOMERS OF Li 10 

We have carried out similar calculations which have led 
us to two geometries of Liio which appear to be serious 
candidates for low- lying isomers. The first of these is 
the very asymmetrical structure shown in Fig. 0] (left). 
Here, with an atom, roughly speaking, at the 'center' 
of the cluster, there can, of course, be no symmetrical 
arrangement around it. 

We have evaluated the total energies, the distance ma- 
trix (in A) and the calculated frequencies (in cm -1 ) of 
the geometrical structure shown in Fig. 0] (left) 5 . The 
second structure we found as a low-lying isomer is shown 
in Fig. Q] (right) 5 . This second isomer is characterized 
by symmetry, and may be thought as obtained from 
the Lis cluster, by 'adding' two more Li atoms (with, 
of course, different bond lengths and angles), shown in 
Fig. Q] (right) in the top- right-hand and bottom- left-hand 
positions. 



V. DISCUSSION AND SOME FUTURE 
DIRECTIONS 

Let us begin this discussion of the results presented 
above by dealing first with further trends, to add to those 
reported in Section Iffl 



A. Sum of eigenvalues (E 3n ) related to total energy 
E n , for n atom cluster 

In Fig. \E\ we have plotted for n — 2 to 10 the total 
energies E n versus the eigenvalue sum E sn from the HF 
occupied eigenvalues ei„, defined by 

E sn = ^ fm- (1) 
occupied i 



FIG. 5: Total energies E n versus eigenvalue sum E 3n of Li n 
(n = 2 - 10). The straight line is Eq. ffi. 

We find an approximately linear relation, which is 
compared with a theoretical prediction from initially 
Thomas-Fermi statistical theory, given by March and 
Plasketl£ for neutral atoms: 

E =\ E s- (2) 

The result @ was found, independently, by Ruedenberg? 
from self-consistent field results for molecules, and its 
foundation in density functional theory was subsequently 
discussed by one of us£. 



B. Trends in dissociation energy D n with number 
of atoms n in Li n clusters denned as \E n — nEi\ 

As a further trend through the clusters Li„ with n go- 
ing from 2 to 10, we have collected in Table [I] results 
from our HF+MP2 calculations for the dissociation en- 
ergy D n . 

Mucci and MarcbS, in early work, stressed the merit of 
Teller's theoremiS, which states that molecules/clusters 
do not bind in any wholly local density approximation 
(LDA). To avoid confusion with current terminology, 
'wholly LDA' in Teller's theorem refers to also treating 
the kinetic energy T by the Thomas-Fermi (TF) result 

r TF = c k / W r)]»/'*, c k = ^(lf\ (3) 

where p(r) denotes the electron density of the 
molecule/cluster. If T denotes the correct (single- 
particle) kinetic energy at the equilibrium geometry, 
Mucci and March— pointed out that the difference be- 
tween T and Ttf was entirely due to electron density 
gradients, e.g. Vp, V 2 p, etc. 
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FIG. 6: Dissociation energy D„/N 2 , where N is the total 
number of electrons in cluster for Li n , versus HF kinetic en- 
ergy jThf' mm us Thomas- Fermi kinetic energy . 

This idea was followed up by Allan et al~, who showed 
that D/N 2 , with TV the total number of electrons in the 
molecule, indeed correlated with the simplest energy con- 
structed from V/9, namely the von Weizsacker (vW) in- 
homogeneity kinetic energy defined by 

r i/2£>!*. (4) 

8 J p 

This is the first-order correction due to density gradients 
in the difference T — Ttf- Therefore, in Fig. we have 
plotted D n /N 2 , already recorded in TableQJ versus T-jjp — 

Tr£p (where = —E^p at equilibrium from the virial 
theorem) for the series Li„, with n again running from 2 
to 10. 

Subsequently, one of usA2> has given the inhomogeneous 
electron gas theory of molecular dissociation energies. 
This reveals that T\{? is the crucial variable to use in 



n 


Symmetry 


D n 


N Dn/N 2 


2 




0.02530 


6 0.000702778 


4 




0.07896 


12 0.000548333 


6 


D 3h 


0.14086 


18 0.000434753 


6 


D 4h 


0.15411 


18 0.000475648 


8 


T d 


0.23765 


24 0.000412587 


8 


c s 


0.24471 


24 0.000424844 


10 




0.31460 


30 0.000349555 


10 


Ci 


0.31843 


30 0.000353811 



TABLE I: Dissociation energy D n = \E n — nE\\, where E n 
is energy of Li„ and E\ that of isolated Li atom. The sec- 
ond column refers to symmetry, when two low-lying isomers 
have been treated, while the fourth column refers to the total 
number of electrons, N. All energies are in Hartree. 



characterizing D n /N 2 . Therefore, in Fig.|§]we have made 
a further plot of D n /N 2 already given in Table QJ versus 

I-^hf — ^tf I' tne ia tter quantity being expected to re- 
flect T^y, albeit approximately. There is a quite clear 
correlation from Fig. between D/N 2 and the 'gradient 
difference' Thf — Ttf, substantiating the work reported 
in Refs. llnEl 

C. Comments on individual clusters and especially 
role of electron correlation 

In cases when n — 6 and 8, two low- lying isomers have 
emerged from our studies using the HF+MP2 approxi- 
mation. A strong pointer to the importance of electron 
correlation in any quite decisive prediction of the true 
ground-state geometry is afforded by the two structures, 
considered for Lis , where Kornath et aZi have experimen- 
tally verified from their Raman studies that the correct 
structure is Td. This is indeed predicted to be lower in 
energy than C s in our HF studies, but the situation is 
(wrongly!) changed over by treating electron correlation 
at MP2 level. However, as noted, HF+MP2 vibrational 
frequencies are considerably higher for T^ than for C s , 
the latter clearly being in marked disagreement with ex- 
periment! On the other hand, a comparison of HF+MP2 
energies and HF energy plus bond-order correlation en- 
ergy fSec lVDjl shows that indeed inclusion of bond-order 
correlation allows to confirm the experimentally observed 
Td structure for Lis. 

Turning more briefly to Lig, we find two low- lying iso- 
mers, both with lower energies than either 6 isolated Li 
atoms or 3 isolated Li2 molecules. Also stability with re- 
spect to isolated L14 and Li2 components is clear. Again, 
however, while lies higher in energy than D 3 ^ at the 
HF level of approximation, including electron correlation 
at order MP2 reverses the ordering. Our conclusion here 
then which seems to us firm is that two low lying isomers 
of Lie are found. However, the energy ordering is not de- 
cisive, though our prejudice here is in favor of the MP 2 
addition, namely D^h symmetry, but that must remain 
conjecture until electron correlation is treated by more 
refined approaches such as the coupled cluster approxi- 
mation. Finally, the largest cluster studied here seems 
fairly 'strongly' bound, being stable again with respect 
to 10 Li atoms, 5 Li 2 dimers, and also the isolated 'frag- 
ments' Lig + Li4 and Lis + Li2- None of this must be 
taken to mean that our structure for Liio is the lowest- 
lying isomer of this cluster, though that is our prediction 
at the HF+MP2 level. 

In the light of the importance of electron correlation, 
we shall now present some 'heuristic' ideas on this sub- 
ject, in which appeal will be made to bond-order ver- 
sus bond-length relations which we have used earlier for 
polyatomic molecule^. We note here that the study of 
Kornath et alX already includes extensive and careful ab 
initio MO calculations with high level of treatment of 
electronic correlations. 
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D. Some approximate considerations on the 
magnitude of electron correlation energy in the Li n 
clusters with n from 2 to 10 

Using the Lowdin dennitioni^*iSii& of electron correla- 
tion, say £7 L , as 

^») = ^L-4f ) ! (5) 

we might use, as a first approximation to the cor- 

relation energy in the K shells, namely 

E^ wn^W. (6) 

One can utilize, for example, the recent study of Alonso 
et alH to estimate Ec for the Li atom ground state. 
Their pairing energy E ss is, for neutral atoms, given in 
their Table 2 as —1.28 eV, and hence for the Liio clus- 
ter we estimate (essentially from 10 separate K shells) 
the correlation energy in magnitude to be 0.47 Hartree. 
The upper curve in Fig. [7J shows this K-shell-like mag- 
nitude of E c for the Li n clusters under consideration. 
Of course, the original 2s electrons will form molecular 
orbitals where again there is pairing of electrons with an- 
tiparallel spins, and we therefore suppose that the upper- 
most curve in Fig. will be below the 'true' magnitude 
of the correlation energy curve. 

For comparison with this K shell estimate, we have 
plotted the 'correlation energy' given by HF+MP2 per- 
turbation theory. The circles in Fig. \7\ show our results 
for Li„ (n = 2,4,6,8 and 10). Two low-lying isomers 
were considered in the present study of both Lig and Lis. 
The small-dashed curve is proportional to n as in Eq. JHJ , 
whereas the long-dashed curve is linear in n but does not 
pass through the origin, in contrast to Eq. JJjJl. 

Evidently, while MP 2 corrections obviously improve 
the HF energies, they yield only ~ 50 to 60 % the Li cor- 
relation energies, the poorest result being for the dimcr. 

Following the technique reported in Ref. Il3l we have 
calculated the correlation energies for lithium clusters us- 
ing the HF bond orders. Fig. [7J compares our results with 
the values obtained from MP2 calculations. In particu- 
lar, the experimental correlation energy i?° xp is defined 



n 




Ef° 


I-.MP2 


17BO 771MP2 


= Ef°/Z 


4 




-0.27126 


-0.09887 


-0.17239 


-0.06782 


6 


D 3h 


-0.41049 


-0.15191 


-0.25858 


-0.06842 


6 


D 4h 


-0.42696 


-0.16724 


-0.25972 


-0.07116 


8 


T d 


-0.57063 


-0.23078 


-0.33985 


-0.07133 


8 


c s 


-0.57561 


-0.24506 


-0.33055 


-0.07195 


10 


D 2d 


-0.72000 


-0.30339 


-0.41661 


-0.07200 


10 


Ci 


-0.72041 


-0.31088 


-0.40953 


-0.07204 



TABLE II: Correlation energies from both bond order (BO) 
and MP2 calculations. The second column lists the symmetry 
of the cluster. All energies are in Hartree. 
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FIG. 7: Different estimates of the magnitude of the electron 
correlation energy in Li n clusters versus n (Sec. ED) . Solid 
line is based on 'estimate' of n 'independent' K shell contri- 
butions and here correlation energy is directly proportional 
to n. Open circles show results due to MP2 perturbation 
theory. Two different 'linear' approximations are represented 
by the dashed lines. Small-dashed line passes through origin 
at n = 0. Long-dashed line is a (slightly) better fit to the 
MP2 results shown as circles. Open triangles show results 
due to bond-order correlation, for which the dash-dotted line 
is a guide for the eye. For Li6, Lis and Liio, values for two 
low-lying isomers studied in each case are shown. 



as: 

£ c cxp =E S - E HF , (7) 

with Es the exact Schrodinger e nerg y and -Ehf the 
Hartree- Fock energy. Following Ref. [jji we calculate the 
theoretical correlation energy Ef° as the sum of atomic 
and bond contributions: 

N JV-1 N 

E*° = E E§n cS + E E a « B *i' ( 8 ) 

i—l i— 1 j=i-\-l 

with N the number of atoms in the cluster. The Eg 
term is the atomic Schrodinger correlation energy that 
can be obtained from experimental data, as explained in 
Ref. ll3t n e ff denotes the effective atomic electron number, 
which takes into account the electrons not involved in the 
bonds, i.e. 



with Z the atomic number and nbond the total number 
of electrons involved in molecular bonds for each atom. 
Bij is the bond order (BO) between atoms i and j, as 
obtained from HF calculations, and is a parameter de- 
pending on the particular bond i—j. For Li clusters, this 
parameter was obtained from the experimental correla- 
tion energy of Li 2 molecule and its value is 0.06157 a.u. 
While up to now all binding and correlation energies have 
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been calculated within HF+MP2, with full geometry op- 
timization, the bond order correlation energies discussed 
here refer to the optimized geometries within HF. 

Table ITTI reports the calculated correlation energies for 
Li clusters obtained from MP2 calculations (E^ 2 ) and 
by using Eq. © (Ef°), along with the HF energies. As 
shown in Table ITTI and in Fig.[7| the calculated correlation 
energies from Eq. © are higher than the MP 2 energies, 
giving for each Li cluster a total molecular energy less 
than that calculated at HF+MP2 level. Moreover, it 
is noteworthy that the difference E^ P2 — Ef° increases 
with increasing number of atoms in the cluster, and that 
the correlation energy per electron, e c — Ef° /Z , is about 
constant for all clusters, indicating that the total corre- 
lation energy is proportional to the total number of elec- 
trons in the molecule. In addition, a comparison of the 
total HF+MP2 and the HF+BO energies in Table ITTT1 
shows that, while for n = 6 both models predict the clus- 
ter with D4/1 symmetry to be stabler than the cluster with 
D^h symmetry, at variance with the HF result, both for 



n = 8 and n = 10 only the HF+BO energies do predict 
the Td and D 2 d to be stabler than the C s and C\, re- 
spectively, in agreement with the purely HF results and, 
in the case of Lis, m agreement with the experimental 
resullpi. 

Fig. [7| shows that the MP2 correlation energy esti- 
mated with the 6-311G* basis set is systematically much 
smaller than the sum of the (Is) 2 atomic pair correla- 
tion. The referee has pointed out to us that this prob- 
ably results from the basic set deficiency for the inner 
shell correlation. 
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n 




-Ehf 


Ebo 


EmP2 






AEhf 
(xlOOO) ( 


ASbo AEmp2 
xlOOO) (xlOOO) 


4 




-29.75973 


-30.03099 


-29.85860 










6 




-44.65841 


-45.06890 


-44.81032 










6 


D 4h 


-44.65633 


-45.08329 


-44.82357 




- D 3h 


2.0800 - 


-14.390 -13.252 


8 


T d 


-59.56615 


-60.13678 


-59.79693 










8 


c s 


-59.55893 


-60.13454 


-59.80399 


c s 


- T d 


7.2191 


2.239 -7.058 


10 


D 2d 


-74.46031 


-75.18031 


-74.76370 










10 


Ci 


-74.45665 


-75.17706 


-74.76753 


Ci 


- D 2d 


3.6587 


3.249 -3.827 



TABLE III: Hartree-Fock and total energies from both bond order (BO) and MP2 calculations, defined as Ebo = Ehf + Ef 
and Emp2 = -Ehf + E^ iF2 . The second column lists the symmetry of the cluster. Last three columns list differences of these 
energies between clusters with same n but different symmetry (D^ — D^h, for n = 6; C s — T d , for n — 8; Ci — T>2d, for 
n — 10). All energies are in Hartree. 



